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ABSTRACT 

We compute the contributions to large-pT J/ip production at HERA coming from 
fragmentation of gluons and charm quarks. We find that the charm quark fragmenta- 
tion contribution dominates over the direct production of J/ijj via photon-gluon fusion 
at laige-pT, while the gluon fragmentation is neghgibly small over the whole range of 
Pt- An experimental study of pt distributions of J/ip at HERA will provide a direct 
probe of the charm quark fragmentation functions. 
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The anomalously large cross-section for J/ip production at large transverse mo- 
mentum, measured [Q] recently in the CDF experiment at the Tevatron has led 
to a revision of earlier ideas based on the lowest order QCD process of parton fu- 
sion. In this approach 0, ^ , the dominant contribution to quarkonium production was 
expected to come from quark- ant iquark or gluon-gluon fusion, leading to the forma- 
tion of a heavy-quark pair in a colour- singlet state with the correct spin, parity and 
charge- conjugation assignments projected out. Several recent works @, H, ^ 0] have 
drawn attention to additional contributions to quarkonium production coming from 
the fragmentation of gluons and heavy quark jets. Even though the fragmentation 
contributions are of higher order in compared to fusion, they are enhanced by pow- 
ers of p^/m? , where m is the heavy quark mass. Consequently, they can overtake the 
fusion contribution at pt ^ Indeed, the CDF J/ip production data |1|] has been 
successfully explained by several authors |^ by taking into account both the fusion 
and fragmentation contributions. The gluon fragmentation contribution is found to 
dominate over fusion at large pt {pt > 5 GeV) , while the charm quark fragmentation 
contribution is much too small. As we shall see in this letter, the photoproduction of 
J ftp at HERA presents a complementary process - i.e. the charm quark fragmenta- 
tion is expected to overtake fusion at large-pr, while the gluon fragmentation remains 
small. Thus a measurement of the large-pr J/ 4^ production cross-section at HERA 
will provide a valuable probe for the charm quark fragmentation contribution. 

A brief discussion of the colour-singlet model, used in the computation of both fu- 
sion and fragmentation contributions, is in order. Strictly speaking, the colour- singlet 
model is a non-relativistic model where the relative velocity between the heavy quarks 
in the bound state is ignored. However, in general, the relative velocity, v, in quarko- 
nium systems is not negligible and 0{v) corrections need to be taken into account. 
Starting from a non-relativistic QCD Lagrangian, a systematic analysis using the fac- 
torisation method has been recently carried out by Bodwin, Braaten and Lepage j^. 
In this formulation, the quarkonium wave-function admits of a systematic expansion 
in powers of v in terms of Fock-space components : for example, the wave-functions 
for the P-state charmonia have the conventional colour- singlet P-state component at 
leading order, but there exist additional contributions at non- leading order in v, which 
involve octet S-state components; i.e. 

\Xj) = 0(l)|QQpPf]) + Oiv)\QQfs[%) + ... (1) 

In spite of the fact that the octet component in the wave-function is suppressed by a 
factor of V, it is important for the decays of P-states for the following reasons : 1) 
The P-state wave-function is already suppressed by a factor of v owing to the angular- 
momentum barrier; but the corresponding colour-octet component is an S'-state which 
is unhindered by this barrier. The colour-octet component can, therefore, easily com- 
pete with the colour-singlet. 2) The second reason is even more compelling : the 
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perturbative analyses of P-wave decays of quarkonia |]TT| reveal a logarithmic infrared 
singularity. But the octet component allows the infrared singularity to be absorbed via 
a wave-function renormalisation, without having to introduce an arbitrary infrared cut- 
off. So a consistent perturbative treatment of P-state decays necessarily involves the 
octet component. The price to pay for this is that two independent matrix-elements, 
viz., the singlet and the octet matrix elements are needed, unlike the case of the 
colour-singlet model where the entire long-distance information could be factorised 
into a single non-perturbative matrix-element. As in the case of the P-state decay 
widths, the P-state fragmentation functions also involve the octet component P| . The 
octet component appears in the computation of the fusion contribution as well, but is 
negligible in the large-p^ region of our interest. 

For S'-state resonances like the J/ip and the ip', the octet contribution is suppressed 
by powers of v. Further, the S'-wave amplitude is not infrared divergent and can, there- 
fore, be described in terms of a single colour- singlet matrix-element. But recently, the 



CDF collaboration has measured [0 the ratio of J/ip^s coming from x decays to 
those produced directly and it turns out that the direct S-state production is much 
larger than the theoretical estimate. It has been suggested |jl3|] that a colour octet 



component in the S'-wave production coming from gluon fragmentation as originally 



proposed in Ref. |T^, can explain this J/ip anomaly. This corresponds to a virtual 
gluon fragmenting into an octet ^Si state which then makes a double El transition 
into a singlet ^Si state. While this process is suppressed by a factor of v"^ as compared 
to the colour-singlet process, it is enhanced by a factor of a^. One can fix the value of 
the colour-octet matrix-element by normalising to the data on direct J/ip production 
cross-section from the CDF experiment. The colour-octet contribution to S'-state pro- 



duction has also been invoked to explain the large ip' cross-section measured by 



CDF [|^, but there can be a large contribution to this cross-section coming from the 
decays of radially excited P-states [|^. Independent tests of the S-state colour octet 
enhancement are important and there have been recent suggestions [0 as to how one 
can use e^e~ collisions to probe the octet contribution. Thus, the possibility of a large 
colour-octet contribution to the S-state fragmentation function remains open, though 
theoretically less compelling than for the P-state. We shall see below that the inclusive 
photoproduction of J/ip is insensitive to the former, but it is sensitive to the latter. 

In this letter, we study inclusive J/ip production in ep collisions at HERA. The 
fusion contribution to the photoproduction of J/ ip in the colour-singlet model [@] comes 
from photon-gluon fusion. Recently, the next-to-leading order corrections to this pro- 
cess have been computed within the colour- singlet model and compared 0] with 



the results on integrated cross-sections from HERA; and it has been found that the 
integrated cross-sections at next-to-leading order are in reasonable agreement with 
the data. The integrated cross-sections are, however, insensitive to the fragmentation 
contributions, because the latter dominate only at large p^- To get a handle on the 
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fragmentation contributions to J/%1) production at HERA it is important to study the 
Pt distributions, rather than integrated cross-sections. 

The fusion contribution to the photoproduction of J/ -0 in the colour-singlet model 
takes place through the following subprocess: 

1 + 9^ ccp5«] + g, (2) 

where the J/ip is taken to be the colour-singlet ^5*1 cc state. The pr differential cross- 
section for the photoproduction of J/ if) in the colour-singlet model is given as 

2 d2 



dcr _ f 128Tr^alapTxG{x)z{l - z)MelR^Q ^ ^ /on 

-J"^' 27\MHl - z) + vl]^ ^'^^ 



where 



dpT J 27[M2(1 -z)+ pI\ 
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(M2+p2)2 [p2^ + M2 (1-^)2; 
4 4 

+ '-Pi (4) 

(M2 + P|)2[p2. + M2(l - Z)2]2 ■ ^ ^ 



In the above equation, the variable z is the inelasticity variable, defined as 

Pj-Pp' 

and X is related to px and 2; , 



(5) 



I r ;\^2 ^2 
X = - 

s 



M2 ^ p2^ 



2; z{l — z) 



(6) 



where s = 4EpU is the photon-proton cm. energy. As usual, M and Rq denote the 
J/ip mass and wave function at the origin. 

The fragmentation contribution is computed by factorising the cross-section for 
the process 7p — * {J/ilj,Xi)X into a part containing the hard-scattering cross-section 
for producing a gluon or a charm quark and a part which specifies the fragmentation 
of the gluon or the charm quark into the required charmonium state, i.e. 

da{-fp {J/ip, Xi)X) j du da{-fp cX)D^_^j/^^^.){uj, /i), (7) 

c "'O 

where c is the fragmenting parton (either a gluon or a charm quark). D(u},iJ,) is the 
fragmentation function and cu is the fraction of the momentum of the parent parton 
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carried by the charmonium stateQ. The fragmentation function is computed pertur- 
batively at an initial scale /iq which is of the order of m^. If the scale /i is chosen to 
be of the order of pt, then large logarithms in fi/rric appear which have then to be 
resummed using the usual Altarelli-Parisi equation: 

o[i ■ Jw y y 

where the Pij are the splitting functions of a parton j into a parton i. We consider the 
fragmentation of gluons and charm quarks alone since the light quark contributions 
are expected to be very small. The gluons are produced via the Compton process: 

l + + (9) 

whereas the charm quarks are produced via the Bethe-Heitler process: 

7 + + (10) 

Using these cross-sections, we compute the fragmentation contribution to da/ dpx which 
is given by a formula similar to Eq. but with an extra integration over cu, or equiv- 
alently over x, because of the relation 

uj = h z. (11) 

xsz 

For the fragmentation functions at the initial scale /i = /ig, we use the results of Refs. 
and P for the gluon fragmentation functions into J/ip or x states, and Refs. and 
for the corresponding fragmentation functions of the charm quark. These fragmen- 
tation functions include the colour-octet component in the P-state, but do not include 
any colour-octet contribution in the S-state. For the case of gluon fragmentation, we 
have separately studied the effect of the S'-state colour-octet component by modifying 



the fragmentation functions as in Ref. |T^. For the charm fragmentation, the S'-state 



colour-octet contributions are sub-dominant and we have neglected these contributions. 
In principle, at HERA energies we can also expect contributions from 5-decays but 



these turn out to be dominant at values of ^ < 0.1 |T9[, and can, therefore, be safely 
neglected in our analysis. 

We have computed the cross-sections for two representative values of the photon 
using the MRSD-' structure functions PO] and we have used Q = pt/^ as 



the choice of scale. In principle, one can integrate over the photon energy spectrum; 
but for the purposes of studying the relative magnitudes of the fusion and fragmenta- 
tion contributions to the cross-sections, it is enough and indeed more transparent to 



^We use the notation u instead of the more usual z to avoid confusion with the inelasticity param- 
eter, defined in Eq. ^ 
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present the results for fixed values of v. In Fig. 1, we present the results for dcr/dpT 
as a function of pt, for z/ = 40 and 100 GeV. For the inelasticity parameter, we use 
the cuts 0.1 < z < 0.9, as used in the ZEUS experiment at HERA We find that 
the fusion contribution, shown by the solid line in Fig. 1, is dominant at low pt, but 
the charm quark fragmentation contribution (shown by the dashed-dotted line) be- 
comes important for values of pr greater than about 10 GeV. The gluon fragmentation 
contribution (shown by the dashed line in the figure) is smaller by over an order of 
magnitude throughout the range of pt considered. Also shown as the dotted line in the 
figure is the gluon fragmentation contribution including the octet contribution for the 
S'-state, where the numerical value of the octet 5'-wave matrix element has been de- 



termined |T^, |14[ from the CDF data ||T2[ as mentioned above. The major uncertainty 
in the prediction is due to the limited information we have on the colour-octet matrix 
elements; the normalisation of the fragmentation contribution can change by a factor 
of 2-3, due to this uncertainty ||^. Moreover, next-to-leading order corrections will 
also change the absolute normalisation of our predictions - for the fusion prediction 
this is expected to give an enhancement (K-factor) upto a factor of 2 [|rS| and similar 
K-factors are also expected in the case of fragmentation contributions. However, our 
choice of Pt/2 as the scale (instead of pt) is expected to account for the i^'-factor 
enhancement, at least in part. 



The charm fragmentation subprocess (Eq. |ToD is gluon-initiated while the gluon 
fragmentation subprocess (Eq. ^ is quark-initiated. This explains why the charm 
fragmentation process dominates. It is important to note that gluon fragmentation 
turns out to be the most important source of J/ip production at the Tevatron, while 
the complementary information on the fragmentation of charm quarks can be studied 
at HERA. An experimental study of pt distributions at HERA will provide us with 
the first direct measurement of the charm quark fragmentation functions. 

Since the J/ip^s produced in the fragmentation process are softer in energy than 
those produced via fusion, it turns out that the average value of z for the former are 
smaller than the latter. To enhance the fragmentation contribution, it is efficient to 
use a stronger upper cut on z. In Fig. 2, we have shown the cross-sections for J/ip^s 
produced via fusion and from charm quark fragmentation with z < 0.5. We find 
that this cut helps to cut down the fusion contribution to J/ip without significantly 
affecting the charm fragmentation contribution, thereby providing a better signal for 
the fragmentation process. 

To summarise, we have studied the pt distribution of J/ ip cross-sections at HERA 
coming from the fusion and fragmentation processes. We find that the large-p^ end is 
dominated by contributions from charm quark fragmentation. The information that 
can be obtained from HERA is thus complementary to that obtained from the large- 
ly J/'ip production at the Tevatron, which is dominated by the gluon fragmentation 
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contribution. By applying judicious cuts on the inelasticity parameter z it is possible 
to enhance the charm fragmentation contribution relative to the fusion contribution. 
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Figure 1: da/dpx (in nb/GeV) for inclusive J/ ip production at HERA for photon energy 
u = 40 GeV (upper figure) and u = 100 GeV (lower figure). The solid line represents 
the fusion contribution and the dashed-dotted line the charm quark fragmentation 
contribution. The dotted and dashed lines represent the gluon fragmentation contri- 
butions with and without a colour-octet component for the (S-state. The cut on the 
inelasticity parameter, z, is 0.1 < z < 0.9. 
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Figure 2: da/dpx (in nb/GeV) for photon energy u = 100 GeV. The sohd hne rep- 
resents the fusion contribution, and the dashed hne the charm quark fragmentation 
contribution, using a cut z < 0.5. 
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